Fluctuations in Y Ba 2 Cu2,OQ^ single crystal: an evidence for 2D — > 3D crossover 
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Magnetization measurements as a function of temperature are reported for Y Ba2Cu2,Oe.5 crystal 
(T c — 45.2 K) for fields between 0.2 to 3.5 Tesla. All isochamps for H > 1 Tesla intersect at 
£NJ ■ T£ D — 42.8 K, implying fluctuations contribution to the magnetization. These curves collapse into a 

' single curve when magnetization and temperature are scaled according to the predicted "2D scaling" 

, in the fluctuation regime. Surprisingly, the low field curves also intersect, at T£ D ~ 43.4 K, and they 

obey a 3D scaling. We provide a theoretical picture of the magnetization in the fluctuation regime 
based on the Lawrence-Doniach model. Within this model we calculate the field and temperature 
Q\ ' dependence of the magnetization. The two intersection points and the 2D — > 3D crossover are 

consistent with the experimental observation. 
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High-temperature superconductors (HTS) are characterized by a wide temperature range in which fluctuations are 
important 0], This range is proportional to the Ginzburg parameter Gi , which is very sensitive both to the 
dimensionality D of the system and to the superconducting coherence length £. Thus, strong fluctuations, usually 
negligible in conventional superconductors, become extremely important due to the small £ and quasi two-dimensional 
structure. 

A useful tool in the analysis of the nature of critical fluctuation is the dimensionality dependent scaling of the 
magnetization M versus temperature T for various DC magnetic fields H H in part of the phase diagram that is "not 
very far" from H c2 {T) Q. The scaled magnetization m = Mj (ht)^ ^ 1 ^ is plotted versus the scaled temperature 
ax — {T — T C (H)) / (HT)^ ^ 1 ^ and all isochamps are predicted to collapse onto a single curve according to a 
dimensionality of the system ^ ^. Once such a scaling is found, the "fluctuation" dimensionality D can be determined. 

y—i \ In 2D systems the relative contribution of the fluctuations to the magnetization is much larger than in 3D systems. 
^ ■ Experimentally, the scaling approach has been widely employed to study the highly anisotropic, quasi 2D, Bi — Sr — 
Ca-Cu-0 (2223) § and (2212) j§ §g], where two-dimensional scaling was shown to work very well. The same 
2D scaling seems to work also for Tl - Ca - Ba - Cu - O (2223) Jk], O, ( 2201) §, @, (2212) j§ Q and (2223) 

^1 jUj, and for Hg - Ba- Ca- Cu- O (1223) @, (1201) @, and (1212)||. 

t—I • The contribution of fluctuations to the magnetization is also borne out in the experiment as a crossing point of all 
isochamps at a temperature T* M, Rj, At this temperature M is independent of H for a large range of fields. This 
feature was previously observed experimentally in 2D systems where 2D scaling was expected [|[ |[ [)[ |To[ [ill |lj. [D^, 
|i~4|| . Theoretically there is evidence that in both 2D and 3D the intersection is not perfect, the intersection points 
however are very close to each other, especially in 2D fl5|| . Though most HTS are either 2D or 3D materials, there 
is, in principle, a possibility that both 3D and 2D behaviors would be measured at the same sample, depending on 
magnetic field and temperature. Such 2D to 3D crossover in vortex fluctuations is expected for highly anisotropic 
superconductors, at high temperatures, simply because the coherence length £ diverges as T approaches the transition 
temperature T c . An evidence for a temperature induced crossover was found in magnetization measurements in 
YDa2Cit307[[l6]]. Another possible experimental approach to study the 2D — > 3D crossover may be based on the 
expected change in the anisotropy caused by changing the oxygen content in Y — Ba — Cu — O. In this system, the 
anisotropy increases with the decrease in the oxygen content. Indeed, a 3D scaling was observed in a fully oxygenated 
YBa2Cu30r single crystal fi~7|| , but a 2D scaling was demonstrated [[l8| in Y Ba^Cu^O^.Q. 

In the present work we establish, experimentally and theoretically, the existence of a 2D to 3D crossover in the 
nature of fluctuations in a high T c superconductor. Specifically, in the oxygen deficient YBa2Cu?,0§.^ (T c « 45.15 K) 
single crystal, at high fields (above 1 Tesla), the magnetization isochamps intersect at one temperature T 2D (Fig.l) 
and exhibit a 2D type of scaling. However, at the lower fields we find another, somewhat smeared, crossing point at 
T£ D (Fig. 2), and the magnetization exhibits a 3D scaling. The presence of two intersection points, as well as 2D and 
3D scaling of magnetization in different field regimes, are the new experimental points in this work. If one defines a 
field dependent "intersection point" as an intersection of two lines for close magnetic fields, one observes that as field 
is lowered, the intersection point first "sits" at (M^jj, T£ D ), than jumps quickly to (M^ D ,T^ D ) and nearly stops there. 
Preliminary discussions of these results were presented in Ref. p9j. In the current paper we provide a theoretical 
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picture of the fluctuation in different regimes of the field-temperature (H — T) plane showing a 2D — > 3D crossover in 
accordance with the experimental results. We calculate the magnetization of the Lawrence - Doniach model describing 
layered superconductors using the " bubble" diagram resummation analogous to that established earlier in the 2D and 
3D limiting cases . The results are compared in the 2D and 3D limits and also between these limits where scaling 
does not hold. 



EXPERIMENT 

Details of sample preparation are given in pl[ . The magnetization of the 2.45 x 3.85 x 0.8 mm 3 Y Ba^GuzO^.^, single 
crystal was measured by a "Quantum Design" SQ UID magnetometer. The high temperature part of the magnetization 
(46 — 200 K) was fitted to a Curie law, M = \H = (xo + G/T)H, and extrapolated to temperatures below T c . The 
extrapolated values of M were subtracted from the raw data measured below T c . This procedure was repeated for each 
value of the applied magnetic field. In Fig. 1 we show the temperature dependence of the magnetization for various 
magnetic fields H > 1 Tesla. All these curves intersect at T = T^ D — 42.8 K, indicating fluctuational contribution 
to the magnetization ||, (?], ||. The subscript 2D is justified by the success of the 2D scaling procedure described in 
Fig. 3 where we plot Mj (HT) versus (T - T c (H)) I (HT) 1 ' 2 for magnetization curves between 1 and 3 Tesla. 

Low-field measurements (H < 1 Tesla) are shown in Fig. 2. Another intersection point, at T = T^ D ~ 43.4 K, is 
found in this field range. This group of curves can be scaled by using the 3D scaling procedure. This is demonstrated in 

2/3 2/3 

Fig. 4 where we plot Mj (HT) versus (T -T C (H)) / (HTy /0 for magnetization curves between 0.2 and 0.75 Tesla. 
These observations imply a 2D — > 3D crossover in the vortex fluctuation regime of our sample. 



THEORY 



In order to find domains of different fluctuation behavior in the H — T phase diagram one has to calculate the 
fluctuational part of the magnetization M defined by the partition function Z : 
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where F is a free energy. In the general case of a layered superconductor with Josephson inter-layer coupling the 
Hamiltonian Hld has the well known Lawrence-Doniach (LD) form J§||/ ; 



Hld = ^A(e F ) / d 2 r 



~iV - —A ) Vv 



+ 2 -V'n+l| 2 + 



(2) 



(t-f)|V^| 2 + ||^| 4 + ^(VxA) 2 
z on 



The magnetic field is assumed to be constant and oriented perpendicularly to the layers (xy) and its fluctuations 
neglected. We use the Landau gauge: A = (0, Hx, 0) . In Eq.(||) ifj n (x,y) is the order parameter in the n-th layer, 
N(cf) is (the 2D) density of states within the layer, ^ a ;, is the in-plane coherence length, (3 - the Ginzburg-Landau 
(GL) coefficient, 7 = (£ c /d) 2 is a dimensionless parameter describing the inter-layer coupling, where d is the inter- 
layer spacing and t = T/T c . This Hamiltonian describes the strong in-plane superconducting fluctuations and their 
inter-layer interactions and can manifest both 3D and 2D behavior in limiting cases as shown below. 

The nonlinear term | -0 rl | 4 in the Hamiltonian becomes very important in the temperature range |1 — t\ ~ Gi at 
the broad region near the transition line, preventing an exact solution of the Hamiltonian in this regime. We are 
able, however, to apply various approximations, as described below, and to show that free energy F can exhibit a 
3D — > 2D crossover as temperature or field are changed. We then are able to approximate the intersection point of 
the magnetization curves in both 2D and 3D limits. 

The magnetic moment of fluctuations described by the Hamiltonian (|^) may be approximately found in three 
limiting cases: (i) The 3D XY model (ii) the 3D lowest Landau levels (LLL) approximation, and (Hi) the 2D 
LLL approximation pp| . Case (i) is not relevant for our experiments since it applies to too low magnetic fields. Cases 
(ii) and (iii) are studied here. In the region of strong fluctuations it is convenient to expand the order parameter in 
terms of the Landau levels eigenf unctions p3 



^(r) = ^>£»a£ 
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C(r) = 
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where N stands for Landau level number and summation index q bears in mind degeneration of the LLL state. 

In the field and temperature ranges of the experiment one can rely on the LLL approximation 0. Therefore, we 
retain in the Hamiltonian only terms with N = : 



q,k 



(5) 



where 



a H 



2eH 
He 



-Cb + t-i 



is a dimcnsionless temperature parameter, and 



a fei ,qi a fe 2 ,92 afe 3 ,93 a &4 ,94 ■ 

^'i9i 



Here 



fcl + fc 2 = fc 3 + fc 4 



?1 + <?2 = 93 + 94 

I{k l ,k 2l h,k A ) = JcPr ^=1,2,3,4^,9, ( r )- 

Since we are interested in vortex liquid phase we will use the renormalized high temperature expansion proposed 
for the 2D and 31? cases by Thouless and coworkers |2(J. The first step is to perform a summation of all the 

"bubble" diagrams. This is equivalent to a kind of mean field approximation in which |^„| 4 — > \ip n \ 2 (l^nl" 2 ^ ■ In 

this approximation the Hamiltonian (0) becomes 



q.k 



a H + 7(1 - cos kd) + - L L 



iV(eF) 



(6) 



q.k 



The dimensionless average 



2 L X L V 



can be found self-consistently by solving the gap equation: 



A = Tdf3 ^ 1 2Tf3eH 1 

2L x LyN(ep) 2-* Off + 7(1 — cos kd) + A JV(e F )7rhc ^( flff + 7 + A) 2 - 7 2 

The magnetization calculated from Eq.(^), using the mean field Hamiltonian (Q), is 

M=- eN ^f»A. (8) 

For convenience, we convert the expression for A into a dimensionless form 

A = g - 5f (9) 
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where 

b _ H _ 2T c /3eH c2 {0) 
flca(O)' 9 N{e F )nhc 

The dimensionless coupling constant g is proportional to the 2D Ginzburg number. 

We consider two limits for which the gap equation can be solved analytically. The first case refers to the domain of 
the H — T phase diagram where 



7 « \ (V( b + t - !) 2 + 46* + (b + t - 1)) (10) 
namely, for the experiments at relatively high fields. In this case, a well known 2D result (see Q) reads: 

A = (Vbi)hn(u); u= b -±^, (11) 

2v bt 



ho{u) = \/ u 2 + 1 - u 

where 

T - TAH) 
b + t-l = 

The magnetization in this limit demonstrates a well pronounced 2D scaling dependence 



V HT V VHT 

The second case refers to the limit 



M ( T-T C {H) \ 

a ho 7= (12) 



j x , 2((6 + *-l) 3 4((fe + *-l) 3 7 1 



namely for the experiments at relatively low fields. In this case, 

A = [(6*) 2/3 ] ( 9 -) ' f 3D (V) (14) 



where 

1/3 



and 



tany= ^ . (15) 



Apparently, the behavior of the magnetization in this case is caused by 31? fluctuations: 



M ( T-TAH) 

« ho TE^m ( 16 ) 



(iJT) 2 / 3 y (HT) 2 / 3 

In both limits one clearly finds a scaling behavior, manifested in Figures 3 and 4. For an intermediate magnetic field, 
however, scaling is not expected even though the LLL approximation is still valid. In this intermediate case the scale 
is provided by the inter-layer spacing d. 
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TABLE I: Parameters used for calculations of theoretical curves in Figs. 1 , 2, and 5. 



Parameter 


Value used 


T c 


45.15 K 


dH c2 1- „ 
iT |T=T C 


4 Tptln IK 


c , (c\\ 


1A 4 


P 

N(e F ) 


n 65^1 
u.uu 9 


d 


5 A 


0BCS 


2.6 x 10 27 erg -2 


N(e F ) 


4 x 10 27 5— y 

erg cm^ 




0.5 A 


m c _ ( g a b \ 2 


850 


m„h \ tc J 





FITS AND DISCUSSION 

The solid lines in Figures 1, 2, and 5 are the theoretical magnetization vs. temperature curves derived from Eq. 
( |iT| ) for the 2D— behavior (Fig. 1), from Eq. ( |l5|) for the 3D behavior (Fig. 2), and from Eq. (||) for the intermediate 
case (Fig. 5). Also, as implied by Eqs. ( p"2| ) and (|l6|), in both the 2D and the 3D cases the magnetization data is 
expected to scale. Both of these formulas obey the respective scaling conditions as demonstarted by the solid lines in 
Figures 3 and 4 derived from Eqs. (|ll|) and (|l4|), respectively. All the "2D" experimental curves intersect at about 
T-2D I T c — 0.948 and the 3D curves approximately intersect at somewhat higher temperature T£ D / T c ~ 0.960. 
Table I summarizes the superconducting parameters used in our analysis to fit the experimental data. 

The parameters were derived in the following way: The transition temperature T c was derived directly from the 
magnetization data. The slope d ^ 2 \t=t c = 4 Tesla/K, yielding H c2 (0) = 180 Tesla, gave the best fit to the experi- 
mental data for both intersection points. Here we use the notation H C 2(0) to denote T c *Zz 2 \t=t c rather than unknown 
upper critical field at zero temperature. The latter is unknown and sometimes is estimated as 70% of this quantity 
as in BCS theory inapplicable to the present case. The coherence length was defined by £ a b(0) = \J <Pq/2t:H C 2(Q). It 
should be noted that here d ^ 2 \t=t c is the mean field theoretical parameter rather than experimentally measured; 
direct measurement of this value is expected to yild a smaller value due to contribution of fluctuations p5| . The value 
of the dimensionless coupling constant g = 0.07, the interlay er- coupling parameter 7 = 0.008, and the magnetization 
M^ D = —2.5 x 10~ 4 emu at the 2D intersection point, determine the ratio N f ep ^ , d, and £ c (0)- Assuming the validity 

of the BCS expression for /3 ~ 8 7 ^|^ one gets a rough estimate of the density of states. Now we discuss the applica- 
bility of the 2D and 3D limits to discribe the regions around the crossing points as was done on Fig.l and 2. The 
inequality [u] for t = T*/T c = .95, so that 1 - 1 » b = .001 simplify into b » 7*(1 -t) = .ITeslaj H c2 {0). Similarly 
the condition of applicability of the 3D limit ( see Eq. (|l3|)) can be simplified into 6<< 1 — t — 7 = hTesla/ H C 2(0). 
The use of the 2D limit in Fig.l is therefore justified for B = l.bTesla or larger, while for B < .75Tesla the use of 
the 3D limit in Fig. 2 is justified. 

To conclude, a crossover between 2D and 3D behavior in the magnetization of Y Ba,2CusO§.§ was observed and 
described theoretically by employing the the Lawrence - Doniach model in lowest Landau level approximation in the 
fluctuation regime. The model yields analytical expressions for two intersection points of the magnetization curves, 
as observed experimentally. One intersection point, for magnetization curves at relatively high fields, is a result of 
fluctuations in the 2D regime. The second intersection point, for relatively low fields, is a result of fluctuations in the 
3D regime. The model also predict scaling of the the magnetization data in the 2D and the 3D regmes, as observed 
experimentally. 
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FIG. 1: Magnetic moment versus temperature in the high-field region (1.5 — 3 Tesla). The solid lines are fits to Eqs. (|]) and 
(|Ll[) obeying the 2D scaling with the parameters listed in Table I. 
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FIG. 2: Magnetic moment versus temperatures in the low-field region (0.2 — 0.75 Tesla). The solid lines are fits to Eqs. (g) 
and (|l5"|) obeying the 3D scaling with the parameters listed in Table I. 
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FIG. 3: Two-dimensional scaling of the high-field data. The solid line is a fit to Eq. (([l 
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FIG. 4: Three-dimensional scaling of the low-field data. The solid line is a fit to Eq. 
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FIG. 5: Magnetic moment versus temperature for various magnetic fields around 1 Tesla. This field range represents an 
intermediate regime where M(T) does not obey neither 3D nor 2D scaling; It is still described, however, in the framework of 
the Lawrence - Doniah model, Eqs. tfy and (g) , as described by the solid lines. 



